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METHOD FOR THE HIGHER-ORDER BLIND DEMODULATION OF A 
LINEAR-WAVEFORM TRANSMITTER 

5 The object of the invention relates to a process for 
the blind demodulation of signals output by several 
transmitters and received by an array made up of at 
least one sensor. 

10 For example, it applies to an array of antennas in an 
electromagnetic context. 

The subject of the invention is in particular the 
demodulation of signals, that is to say the extraction 
15 of the symbols {a k } transmitted by a linearly modulated 
transmitter . 

Figure 1 shows an antenna processing system comprising 
several transmitters Ei and an antenna processing 
20 system T comprising several antennas R± receiving from 
radio sources at different angles of incidence. The 
angles of incidence of the sources or transmitters may 
be parameterized either in ID with the azimuth 0 m , or in 
2D with the azimuth angle 6 m and the elevation angle A m . 

25 

Figure 3 shows schematically a modulation/demodulation 
principle for the symbols {a k } output by a transmitter. 
The signal propagates via a multipath channel. The 

30 transmitter outputs the symbol a k at the instant k.T, 
where T is the symbol period. The demodulation consists 
in estimating and in detecting the symbols in order to 
obtain, at the output of the demodulator, the estimated 
symbols a k . In this figure, the train of symbols {a k } is 

35 linearly filtered upon transmission by a transmission 
filter H, also called a wave-shaping filter h 0 (t). 

In the rest of the description, the expression "blind 
demodulation" is understood to mean techniques that 



basically use no information on the signal output, 
examples being a wave-shaping filter, a learning 
sequence, etc. 

5 The last ten years have seen the development of SIMO 
(single-input, multiple-output) blind demodulation 
techniques called subspace techniques using 2nd-order 
statistics, as described in reference [7] . However, 
these algorithms have the drawback of not being robust 

10 to either underestimation or overestimation of the 
order of the propagation channel, resulting in temporal 
spreading dependent on the multipaths and on the wave- 
shaping filter- A linear prediction technique, 
described in reference [11], has been proposed for 

15 overcoming this problem, but this has the drawback of 
being less effective when the length of the channel is 
known. To improve the subspace techniques, the method 
described in [16] proposes a parametric technique, but 
unfortunately this requires knowledge of the wave- 

20 shaping filter. 

In reference [13], the authors propose a technique 
based on covariance matching, but this has in 
particular the drawback of being very difficult to 

25 implement. An easier but suboptimal technique, 
described in reference [12], was therefore developed by 
minimizing a likelihood criterion and assuming the 
symbols to be Gaussian in character. This assumption is 
not verified for the widely used linear modulations 

30 such as PSK (Phase Shift Keying) or QAM (Quadrature 
Amplitude Modulation) . 

It is also known, in CMA (Constant Modulus Algorithm) 
methods, to use a spatio-temporal approach described 
35 for example in reference [6] . However, this family of 
methods has the drawback of being suitable only for one 
particular class of modulations, such as PSK, which are 
constant-modulus modulations. This method is iterative 
and therefore has the drawback of having to be 



correctly initialized. Finally, the CMA methods have 
the disadvantage of converging more slowly than the 
abovementioned subspace method. Moreover, reference 
[20] describes a subspace method making use of higher- 
5 order statistics for non-minimum-phase FIR (finite 
impulse response) channel identification. 

The subject of the present invention is a process based 
in particular on blind source separation techniques 

10 known to those skilled in the art and described for 
example in references [4], [5], [15] and [19] assuming 
that the symbols transmitted are statistically 
independent. To do this, the process constructs a 
spatio-temporal observation whose mixed sources are 

15 symbol trains from the transmitter. Each symbol train 
is for example the same symbol train but shifted by an 
integral number of symbol periods T. 

The invention relates to a process for the blind 
20 demodulation of a linear-waveform source or transmitter 
in a system comprising one or more sources and an array 
of sensors and a propagation channel, said process 
being characterized in that it comprises at least the 
following steps: 
25 • the symbol period T is determined and samples are 

taken at T e such that T = IT e (I being an integer); 

• a spatio-temporal observation z(t), the mixed 
sources of which are symbol trains from the 
transmitter, is constructed from the observations 

30 x(kT e ); 

• an ICA-type method is applied to the observation 
vector z(t) in order to estimate the L c symbol trains 
{a m _i} that are associated with the channel vectors 
fi Z/j =fi z (kj) ; 

35 • the L c outputs (a m ,j, fi 2 ,j) are arranged in the same 

order as the inputs (a m -i, h z (i)) so as to obtain the 
propagation channel vectors fi 2 , j=fi z (kj) ; and 

• the phase a imax associated with the outputs is 
determined. 
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The process according to the invention offers in 
particular the following advantages: 

• it makes no assumption about the symbol 
constellations, unlike the methods described in the 
prior art; 

• it requires no knowledge of the wave-shaping 
filter; 

• the modulus of the symbols is not assumed to be 
constant; 

• it is robust to channel length overestimation; 

• it can handle propagation channels with correlated 
paths; and 

• it is direct and simple to implement, with no 
correlated-path crosscheck step. 

Other features and advantages of the subject of the 
present invention will become more clearly apparent on 
reading the following description, given by way of 
20 illustration but implying no limitation, and on 
examining the appended figures, which show: 

• figure 1, an example of an architecture; 

• figure 2, the angles of incidence of the sources; 

• figure 3, the linear modulation and demodulation 
25 process for a symbol train; 

• figure 4, the diagram of a linear-modulation 
transmitter; 

• figure 5, a summary of the general principle 
employed in the invention; 

30 • figure 6, the representation of a constellation; 

• figure 7, a first example of the implementation of 
the method, in which the signal is received baseband; 

• figure 8, a second example in which the signal is 
received in baseband and the multipaths are 

35 decorrelated; and 

• figure 9, a third example in which the signal is 
received in baseband and the multipaths are groupwise 
decorrelated. 
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To explain the process according to the invention more 
clearly, the following description relates to a process 
for the higher-order blind demodulation of a linear- 
waveform transmitter in an array having a structure as 
described in figure 1 for example. 

Before explaining the steps for implementing the 
process, the model of the signal used will . be 
described. 

Model of the signal output by a source or transmitter 
Linear modulation 

Figures 3 and 4 show the process for the linear 
modulation of a symbol train {a k } at the rate T by a 
wave-shaping filter h 0 (t). 

The comb of symbols c(t) is firstly filtered by the 
wave-shaping filter ho(t) and then transposed to the 
carrier frequency fo. The NRZ filter, which is a time 
window of length T, very often defined by 
ho(t) = n T (t-T/2), is one particular nonlimiting example 
of a transmission filter. In radio communication, it is 
also possible to use a Nyquist filter, the Fourier 
transform of which, h 0 (f) * n B (f-B/2), approaches a band 
window B when the roll-off is zero, therefore 
h 0 (f) = n B (f-B/2) (the roll-off defines the slope of the 
filter away from the band B) . 

The modulates signal s 0 (t), output by the transmitter, 
may be expressed at time t k = kT e (T e being the sampling 
period) as a function of the comb of symbols c(t): 



Let the symbol time T be equal to an integer number of 
times the sampling period, i.e. T = IT e , and let 
k = ml+j where 0 ^ j < I. Since c(t) = S r a r 5 ( t-rIT e ) , in 
other words c(t) = a u for t = u!T e and c(t) = 0 for t ^ 




(1) 
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uIT e/ the only values for i for which c((k-i)T e ) is non 
zero satisfy the equation k-i = ul, that is to say such 
that i = ml+j-ul = nl+j , where n = m-u. Finally, 
equation (1) becomes: 

So{mlT e +jT e )= £ h 0 (nlT e +jT e ) am-n for 0<j<l (2) 

5 

The parameter L 0 is the half-length of the transmission 
filter which is spread over a duration of (2L 0 +l)IT e . In 
the particular case of an NRZ transmission filter, 
L 0 = 0. As regards the transmitted signal s(t), this 
10 satisfies the equation s(t) = s 0 ( t ) exp ( j 2nf 0 t ) as it is 
equal to the signal s 0 (t) transposed to the frequency 
f 0 . Under these conditions, the expression s(mIT e +jT e ) 
is, from (2) : 

s(mlTe+jT e )= £ ho(nlTe+jT e )expG27ufo(nl+j)Te)a m .nexpG27cfo(m-n)ITe) 

= £ h F0 (n I T e + jT e ) b m ^ such that 0 < j < I (3) 
where hFo(n" e )=h(iT e ) exp(j27cf 0 iT e ) and bj=aj exp(j27ifoilT e ) 



Reception of the signals by the sensors 

The transmitted signal s(t) (figure 3) passes through a 
20 propagation channel before being received on an array 
made up of N antennas. The propagation channel may be 
modelled by P multipaths of angle of incidence 0 P , delay 
T p and amplitude p p (1 ^ p ^ P) . At the output of the 
antennas is the vector x(t), which corresponds to the 
25 sum of a linear mixture of P multipaths and noise, 
assumed to be white and Gaussian. This vector of 
dimension N * 1 is given by the following expression: 

x(t) = X Pp a(9 p ) s(t-x p )+ b(t)=A s(t) + b(t) (4) 

where p p is the amplitude of the pth path, b(t) is the 
30 noise vector, assumed to be Gaussian, a(0) is the 
response of the array of sensors to a source with angle 



20 
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of incidence 9, and A = [a (0i) . . .a (0 P ) ] and s(t) = [s(t- 
Ti) . . .s(t-tp) ] T . Noting that x p = r p T+Ai p (where 0 ^ Ax p < T 
= IT e and r p is an integer), and inserting equation (3) 
into equation (4), the vector received by the antennas 
5 is given by: 

x(mlT e +jT e ) = £ £ Pp a ( e p) h Fo(n ITe+jTe-Atp) b + b(mlTe+jT e ) (5) 

By making the change of variable according to u p = n+r p , 
the vector received by the antennas is given by: 

x(mlT e +jT e ) = £ I Pp a(9 p ) h F0 ((u p - r p )IT e +jT e -ATp) b m ^ + b(mlT e +jT e ) (6) 

10 Notating r m i n = min{r p } and r max = max{r p }, equation (6) 
may be rewritten as follows: 

X(mlT e +jT e ) =2, Z PP a ( G p) M(U- r p) IT e + jTe-Ax p ) Indj^LO. rp +L 0](U) b m _„ 

+ b(mlT e +jT e ) (7) 

where Ind[ r ,q] (u) is the usual indicatrix function 
(Ind[ r ,q] (u) = 1 for r ^ u ^ p and Ind [r , q) (u) = 0 
15 otherwise) defined over the set of integers relating to 
the value in the binary set {0,1}, characterized by 
Ind[ r ,q] (u) = 1 if u belongs to the interval [r,q] and 
Ind [r ,q] (u) = 0 otherwise. Thus, denoting the channel 
vector by v(t) : 

/' 

v(ulT e +jT e ) = £ p p a(9 p ) h F o((u- r p )IT e +jT e -At p ) lnd Irp -Lo, rp + Loj(u) (8) 



where t = uIT e +jT e and equation (5) becomes: 

x(mlT e +jT e )= X v(ulT e +jTe)b m . u + b(mlTe+jTe) (9) 



Inter-symbol interference 

25 

The observation vector x(t) coming from the antenna 
array at t = mIT e +jT e involves, from equation (9), the 
symbol b m , but also the symbols b m - u where u is a 
relative integer lying within the [r m i n -Lo, r max +L 0 ] 
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interval, which phenomenon is more widely known as ISI 
(inter-symbol interference) . Let L c be the number of 
symbols participating in the ISI and let the interval 
of values taken by the latter be limited. From equation 
5 (9), if the intersection of the intervals [r p -L 0 , r p +Lo] 
is non-empty, then L c = Ir^-r^J + 2L 0 + 1- 
Consequently, when r max = r^n, that is to say when all 
the multipaths are correlated, the lower bound of L c is 
reached and is given by L c = 2L 0 + 1- This case is also 

10 mathematically expressed by | max {x } - min {x } |<T. On 

p p 

the other hand, if the intersection of said intervals 
is empty and if all the intervals [r p -L 0 , r p +L 0 ] are 
disjoint, then L c = Px(2Lo+l), which constitutes an 
upper bound for the set of values that can be taken by 

15 L c . The latter situation corresponds specifically to 
the case of multipaths that are all pairwise 
decorrelated. This may also be mathematically expressed 
as Vi # j,| n - rj l>2L 0 , this condition being obtained 
whenever I t± - Xj | > (2L 0 + 1)T. To summarize, the 

20 quantity L c is in general bounded as follows: 

2L 0 +1 <U<Px(2L 0 +1) (10) 

The equation expressing the vector received by the 

sensors can then be rewritten in the following manner, 

25 but this time only the L c symbols b m _ u of interest 
appear: 

X(mlT e +jTe)=X h(n(l)IT e +jT e )b m .n(i) + b(mIT e +jTe) (11) 

/=l 

where V 1 ^ I ^ L c and r^n-I^o < n(I) < r m i n +Lo and where: 

h(t)=t P P a(9p)h F0 (t-Tp) (12) 

h 

30 

ICA techniques 

The process uses ICA techniques based on the following 
model, given by way of entirely nonlimiting 
illustration: 
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L 

u k =X 9i SiK +n k =Gs k +n k (13) 

where Uk is a vector of dimension Mxl received at time 
k, s ik is the ith component of the signal s k at time k, 
n is the noise vector and G = [gi-.-giJ- The objective 
5 of the ICA methods is to extract the I = L components 
Sik and to identify their signatures gi (the vectorial 
response of source i through the observation u k ) on the 
basis of the observation u k . The number I =L of 
components must not exceed the dimension M of the 

10 observation vector. The methods of references [4], [5] 
and [15] use 2nd- and 4th-order statistics of the 
observations u k . The first step uses 2nd-order 
statistics for the observations u k (these observations 
may be functions of the signals received by the 

15 sensors) in order to obtain a new observation z k , such 

that : 

L 

z k = W 1 Uk=2] gi s ik + n k =6s k + n k (14) 

where the signatures g s (l^i^L) are orthogonal, 
G =[gi-..g L ] and s k = [si k . - - s Lk ] T . The second step 
20 consists in identifying the orthogonal base of the G 
values using 4th order statistics of the whitened 
observations z k . In this way, the signals s k may be 
extracted by effecting: 

Sk = 6 # ZK=6 # W lU K (15) 

25 where s k is the estimate of the signals s k and where # 
is the pseudo-inversion operator defined by 
G # =(G H G)" , G H . 

The ICAR method [19] uses only 4th-order statistics to 
identify the matrix G = [q 1 ... g k ] of signatures. 

30 

To summarize, the idea employed in the process 
according to the invention is to construct a spatio- 
temporal observation, the mixed sources of which are 
symbol trains from the transmitter. Each symbol train 



- 10 - 

is for example the same symbol train but shifted by an 
integral number of symbol periods T. 



Several ways of implementing the method will be 
5 described below, some of which are explained by way of 
non-limiting illustration . 

First way of implementing the process. 

10 Figure 7 shows a first illustrative way of implementing 
the process in which the signal is received in 
baseband. 

The method comprises a step 1.1 of determining the 
15 symbol time T e , for example by applying a cyclic 
detection algorithm, such as that described for example 
in [1] [10] . 

The next step 1.2 consists in interpolating the 
20 observations x(t) with I samples per symbol, such that 
T = IT e . 

Under these conditions where f 0 = 0 and b k = a k , 
equation (11) for the vector becomes: 



Since equation (16) is valid for 0 ^ j < I, the method 




h(n(l) I V jT e ) anwm + b(ml V jT e ) for 0 < j < I (16) 



constructs the next spatio-temporal 
(step 13) from the observations x(kT e ): 



observation 




with h nJ = h(nl T e + ]T e ) and b*(ml T e )=[ b(ml T e ) T ... b(ml T e +( M)T e ) T ] T 



Since it is known that x(t) has the dimensions Nxl, the 
vector z(t) has the dimensions NIxl. 
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h(k) is a vector whose nth component is the kth 
component of the filter that linearly filters the 
symbol train {a m } on the nth sensor. The filter for the 
vector coefficient h(k) depends both on the wave- 
5 shaping filter and on the propagation channel. 

To extract the L c symbol trains {a m -i} of interest (the 
number of symbols participating in the ISI), the method 
samples the received signal with I = (2L 0 +1) , assuming 
10 that P < N. 

Since it is known that the NRZ filter satisfies 2Lo+l = 
1 and the Nyquist filter 2L 0 +1 = 3 for a roll-off of 
0.25, the symbol trains may be extracted for these two 
15 wave-shaping filters when P ^ NI and 3P ^ NI, 
respectively. 

Having determined the observation vector z (t) , the 
process applies an ICA-type method to estimate the L c 
20 symbol trains {a m -i} associated with the channel vectors 

The jth output of the ICA method gives the symbol train 
{^mj} associated with the channel vector . The 

25 estimated symbol trains {a mj } arrive in a different 
order from that of the {a m -i} trains satisfying: 

a mij =pexp(jai)am-i and h zy =h 2 (i) (18) 

The symbol trains {a mj } are estimated with the same 

30 amplitude because the symbol trains {a m _i} all have the 
same power, satisfying the equation: 

E [ I a m -n(i) I 2 ] = - . • =E [ | a m _n<Lc) I 2 ] - 

The next step 1.4 of the process has the objective of 
35 ordering the L c outputs (a mj h^j) in the same order as 

the inputs (a m _i, h z (i)) so as to obtain the channel 
vectors h^j=h z (kj). To do this, the process 
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intercorrelates pairwise the outputs and a^j , 

calculating the following criterion Cj.,j(k): 

., (19) 

When the function |Cj.,j(k)| is a maximum in k = k max , the 
5 ith and jth outputs satisfy the equation: a^ = a m _ kmaxJ . 

The algorithm for classifying the outputs a mn(1) ...a mn(Lc) is 

for example composed of the following steps: 

Step A.l: Determination of the output a mimax associated 

with the channel vector of higher-modulus h z jmax • 
10 Step A. 2: For all the outputs a m _ kj where j ¥= i max/ 

determination of the indices k = kj maximizing the 
I Ci max ,j {k) | criterion. From this is deduced, for each j, 
that a mimax = a m _ kjj . Since it is known that c im ax,j(kj) = 

exp ( jcximax" j j ) the jth output is reset to the same phase 
15 as the iith imax output by taking a m _ kj = c imax , j (kj) a m } . 

The channel vectors are also reset in terms of phase by 

A 

taking: h z (kj) = h zj c imax jCkj)*. 

Step A. 3: This step reorders the outputs a m _ kj and the 

A 

A A 

channel vectors h z j=h z (kj) in the increasing order of 
2 0 the Kj, since it is known that a m = a mimax and that 

After these three steps, the symbol trains {a m - k } 
associated with the channel vectors h 2 (kj) are obtained. 

25 Since it is known that the estimated symbols satisfy 
the equation a ra - k = exp ( ja imax ) a m - k/ the last step of the 
process consists in estimating this phase a imax . To do 
this, the constellation of symbols a k is firstly 
identified among a database consisting of the set of 

30 possible constellations. This database consists of 
known constellations such as nPSK, n-QAM. Each time 
that a new constellation is detected or becomes known, 
this is added to the database. 

Figure 6 shows an example of an 8 -QAM 

35 constellation when aimax=0 and aimax ^0 . In this 



implementation example, the process then includes the 
following steps: 

The next step 1.5 consists in determining the 
output phase associated with the channel vector of 
5 higher modulus. To identify the constellation and 
determine the phase, the process performs, for example, 
the following steps: 
Step 1.5 = Steps B.l, B.2 and B.3 

Step°B.l: Estimation of the positions of the states of 

10 the constellation (red points in the figure) by seeking 

the maximum of the 2D histogramme of the points 

Mk= (real (ak) , imag(ak))- For a constellation consisting 
of M states, M pairs (u m/ v ffl ) for 1 < m < M are 

obtained. 

15 Step°B.2 : Determination of the type of constellation 
by comparing the position of the states (u m /V m ) of the 

constellation of {a*} symbols with a database 
comprising the set of possible constellations. The 
closest constellation is made up of the states (u m ,v m ) 

20 for 1 < m < M. 

Step B.3: Determination of the phase a imax by minimizing 
in the sense of the least squares the following system 
of equations : 

u m =cos (a imax )u m - sin(a imax ) v m 

25 and v m = sin(a imax )u m + cos(a imax )v m for 1 < m < M. 

The process may include a step of estimating the 
propagation channel parameters of angle 0 P and delay x Pr 
of equation (8) by the algorithm proposed in [8] . The 
30 step consists in extracting firstly the vectors h(nIT e + 
jT e ) for 0 < j < I from the channel vectors h r (w y .) 

defined in equation (17). Followed by construction of 
the matrix H= [h (n (1) IT e ) ... h (n (L c ) IT e ) ] from equation 

(11) with the h(nIT e + jT e ) values in order to apply the 
35 parametric estimation method [8] for the multipaths: 

(6 P , x p ) 1 < p < P. 
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Second way of implementing the process 

Figures 8 and 9 show schematically another way of 
implementing the process, which may include two 
variants corresponding to the decorrelated multipath 
5 case and to the groupwise correlated multipath case, 
respectively. 



Decorrelated multipath case. 



10 The signal is received in baseband with {b k }={a k }. 



The multipaths, the delays of which satisfy the 
relationship lij-Xi I > (2L 0 +1)T, have the advantage of 
being decorrelated with one another, satisfying 

15 the equation : E [s (t-Tj.) s (t-Xj) *] = 0. By examining 
equation (4), it may therefore be seen that it is 
sufficient to apply an ICA type method when P < N to 
the observation x(t) in order to obtain the signals 
s(t-ip) for each of the multibars . After estimating the 

20 signals for the various multipaths, the process 

determines their powers in order to keep the signal 

s(t-x pmax ) of the multipath of higher amplitude p pma x. 

This main path is determined using the fact that the 

outputs of the ICA methods asymptotically satisfy: 

p p 
x(t) = £ P P a(0p)s(t-T p ) = £ 315.(0 with (20) 
p = i /=i 

S(t-T ) . 

*,(0 = r=— and ai =yjY p P P a(9 p ) 

y/y p 

25 where y p =p p 2 E [ | s ( t-x p ) 1 2 ] . Since the vectors a(9 p ) are 
normed, satisfying the equation a (G p ) H a (9 P ) =N, the path 
of maximum amplitude will be associated with the i ma x th 
output where ai max =a ini a X H ai m ax is a maximum. From equation 
(3), the output 5 imax (0=s (t-T pmax ) satisfies the 

30 equation : 



^imax 



(mIT e +jT e )= £ h F0 (nITe+jT e -T pmax )a m -n (21) 
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such that 0 < j < I 



10 



and it is possible to constitute the following 
observation vector: 



z (mIT e ) = 



S ims »(mIT e + T e ) 
J iaa ,(mIT e ^(I-l)T e ) 



= X h z (n) a m _ n 



n = - L 9 



(22) 



where h z (n) = 



Vi 



l n,I-\ 



and 



where h n , j=h F0 (nIT e + jT e -t pmax ) . According to the model of 
equation (22), it is sufficient to apply an ICA method 
to the observation z (mIT e ) in order to estimate the 
2L 0 +1 symbol trains {a m - n } with -L 0 < n < L 0 - To extract 
the angles of incidence 8 P of the propagation channel, 
it is sufficient from equation (20) to find, for each 
signature ai (1 < i < P), the maximum of criterion 
c(9) = |a(0) H ail 2 . To extract the delays ii -Xi of the 
propagation channel, it is sufficient from equation 
(20) to find, for each signal (1 < i < P) , the 

maximum of the c (x) = | s^t-r) ^(/)* I 2 criterion. 



15 To summarize, this variant comprises, for example, the 
following steps: 

Step II. a. 1: Determination of the symbol period T, 
applying a cyclic detection algorithm as in [1] [10] . 
Step II. a. 2: Sampling of the observations x(t) with I 
20 samples per symbol such that T = I T e . 

Step II. a. 3: Application of an ICA method to the 
observations x(t) in order to obtain $,(/) and a ± for 

1 < i < P. 

Step II. a. 4: Determination of the output i = imax where 
25 <Xi=ai H ai is its maximum. 

Step II. a. 5: Formation of the observation vector z(t) 
of equation (22) from the signal s imax (t) . 

Step II. a. 6: Application of an ICA method for 
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estimating the symbol trains {a m - n } where -L 0 < n < L 0 . 
From the symbol trains is chosen that one which is 
associated with the higher-modulus vector h z (n), namely 
{a m } . 

Step II. a. 7: Determination of the phase a imax of the 
output associated with the higher-modulus vector h z (n) 
applying steps B.l, B.2 et B.3. 

Step II. a. 8: Phase-resetting of the symbol train {a m } 
by taking a m = a m exp (- ja imax ) - The symbol train {a m } 
constitutes the output of the demodulator of this 
subprocess . 

Step II. a. 9: Estimation of the propagation channel 
parameters, namely angle 9 P and delay x p , by maximizing, 
for l<i<P the |a(9) H ail 2 and | s,(f - r) s x (t) * | 2 

criteria for the angles and delays respectively. 

General case for any or groupwise-correlated multipaths 

In this variant, the diagram for which is given in 
figure 9, the process considers that some of the 
multipaths are correlated. Considering that the 
transmitter is received according to Q groups of 
correlated multipaths, the signal vector received by 
the sensors becomes, from equation (4) : 



Where A q = [a (Gi, q ) ... a(G Pq , q )], Q q =diag ( [p 1#q ...p Pq , q ] ) and 

S (t,X q ) = [S (t-Ti /q )... S (t-Tp q , q ) ] T with Iq=[Ti, q ...T Pq , q ] T . The 

following signals and signatures are estimated as 
output of the separator by applying an ICA method: 



x(t) 




(23) 
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A=[ai . . . apQ,Q]=[AiUi . . . AqUq] n and s(/) =n 



"V, s(M,)" 






V e s(M e ) 







(24) 



10 



where n is a permutation matrix, U q V q - Q q and 



V q E[s (t,T q ) s (t,x q ) H ] V q H =I Pq . Thus, the paths decorrelated 

such that E [s (t~T P/q ) s (t-Tp' /q O *] =0 are received on 
different channels s L (t) and s^t) . The correlated 

paths where E [s (t-x p , q ) s (t-T P ', q ' ) *]*0 are mixed in the 
same channel s.(t) and are present on P Q at the same 

time. In the 1st step of this subprocess, we use this 
result to identify the Q group of correlated 
multipaths. Taking the outputs i and j of the 
separator, the two following hypotheses may be tested: 



Ho 



[5,(0=6,(0 

[Sj(t)=bj(t) 



and Hi 



(25) 



where E[bi(t) bj(t-x)*] = 0 whatever the value of t. 
Thus for the Ho hypothesis, no multipaths exist common 
15 to the two output i and j, and for the Hi hypothesis 
there is at least one of them. The test consists in 
determining whether the outputs s i ,(t) and Sj(t-T) are 

correlated for at least one of the t values satisfying 
Ul < x max . To do this, the Gardner test [3] is applied, 
20 which compares the following likelihood ratio with a 
threshold: 

v^Kind-^L) WW, ^jpme-ty < 26 > 

or Vij(T)<n => hypothesis H 0 
And Vy(T)> n => hypothesis Hi 

The threshold r| is determined in [3] in relation to a 
25 chi square law with 2 degrees of freedom. The output 



associated with the 1st output are firstly sought by 
starting the test by 2 < j < P Q x Q and i = 1. Next, 
removed from the list of outputs are all those 
associated with the 1st which will constitute the 1st 
5 group with q = 1. The same series of tests is then 
restarted with the other outputs not correlated with 
the 1st output in order to constitute the 2nd group. 
This operation will be carried as far as the last group 
where, in the end, no output channel will remain. After 
10 the sorting, what will finally be obtained are: 

Aq=AqU q and s q (/) = V q s(t, i q )for (1<q<Q) (27) 



15 



The angles of incidence 9 p , q are determined from the A q 
values for (1 < q < Q) applying the MUSIC [1] algorithm 
to the A q A q H matrix. The matrices A q are deduced from 
these goniometry values. Since is it known that x q (t) = 
A q s q (t) = A q Q q s(t, x q ) , s(t,x q ) is deduced therefrom to 

within a diagonal matrix by taking s (t,x q ) = A q # x q (t) . 
Since the elements of s (t,x q ) are composed of the 
signals s (t-x p , q ) , the delays x p , q -Ti,i are determined by 
maximizing the c(x) = I £ q , p (t-x) s ±fl (t) *| 2 criteria where 
i q , p (t) is the p th component of s (t,x q ) . 

Since it is known that E[s q (t) s q (t) H ] = I Pq , that A q H A q = 
N I Pq and that A q s q (t) = A q Q q s(t,T q ), it is deduced 
therefrom that the group of multipaths associated with 
25 the largest amplitudes Q q maximizes the following 
criterion: cri(q) = trace (A q H A q ) . From this is deduced 
the best output associated with Aqmax and iqmax(t). Since 



20 



from equation 
equation: 



(3) the vector s(t, 



Zqmax ) 



satisfies 



30 



S(mlT e +jT e ,Iqmax)= 



s(mK e +jT e 
s(mIT e +jX-T qmaP ^) 



= J h F o(nlTe+jTe,Iqmax) a m . n (28) 



n = -L Q 



for 
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0<j<l and with h F o(nlTe+jT e> _lqmax)= 



h F0 (nrT e +fT e -T qmiix] ) 



it is possible to constitute the following observation 
vector from equation (27) : 



z(mlT e )= 



= £ h z (n) a m . n where h 2 (n)= 



VO 



V/-1 



(29) 



5 where h n# j = Vq^* h F0 (nIT e + jT e/ lqmax) - From the model of 
equation (29), it is sufficient to apply an ICA method 
to the observation z (mIT e ) in order to estimate the 2L 0 
+ 1 symbol trains {a m - n } such that -L 0 < n < L 0 . 

10 To summarize, this variant comprises the following 
steps : 

Step Il.b.l: Determination of the symbol period T by 
applying a cyclic detection algorithm as in [1] and 
[10] . 

15 Step II. b. 2: Sampling of the observations x(t) with I 
samples per symbol such that T=IT e . 

Step II. b. 3: Application of an ICA method to the 
observations x(t) in order to obtain s(t) and A from 

equation (24) . 

20 Step II. b. 4: Sorting of the outputs according to Q 

groups of correlated multipaths in order to obtain A q 
and s q (t) for (1 < q < Q) : to do this, a correlation 

test for all the output pairs i and j with the two- 
hypothesis test of equation (26) . Firstly the outputs 
25 associated with the 1st output will be sought by 
starting the test for 2 < j < P Q * Q and i = 1. Next, 
removed from the list of outputs are all those 
associated with the 1st that will constitute the 1st 
group with q = 1. The same series of tests is repeated 



with the other outputs that are not correlated with the 
1st output in order to constitute the 2nd group. This 
operation is continued to the last group where in the 
end no output channel will remain. 
5 Step II. b. 5: Determination of the better group of 
multipaths where (A q H Aq) is a maximum in q = qmax. 
Step II. b. 6: Constitution of the observation vector 
z(t) of equation (29) from the signal s qmax (t) . 

Step II. b. 7: Application of an ICA method for 
10 estimating the symbol trains {a m _ n } where -L 0 ^ n < L 0 . 
From the symbol trains is chosen that one which is 
associated with the higher-modulus vector h z (i), namely 
{ a m -i } . 

Step II. b. 8: Determination of the phase a imax of the 
15 output associated with the higher-modulus vector h 2 (i) 
applying steps B.l, B.2 and B.3. 

Step II. b. 9: Phase-resetting of the symbol trains {a m } 
by taking a m = a m exp (- ja imax ) . The symbol train { a m } 
constitutes the output of the demodulator of this 

20 subprocess. 

Step II.b.10: Estimation of the propagation channel 
parameters, namely the angle 6 q , p and the delay T q , p . The 
angles of incidence 9 q , p are determined from the A q 
values for (1 < q < Q) applying the MUSIC [1] algorithm 

2 5 to the matrix A q A q H . The matrices A q are deduced from 
these goniometry values in order to deduce therefrom an 
estimate of s(t, x q ) taking s (t, x q ) = A q # x q (t). Since the 
elements of the s (t,x q ) are composed of the signals 
s(t-T q , p ), the delays i p , q -Ti,i are determined by 

30 maximizing the c (t ) =|s q , p (t-x ) Si, x (t) *| 2 critera where 
s q , p (t) is the p th component of s (t,T_q) • 

Another way of implementing the process 

Estimation of the carrier frequency and deduction of 
35 the {a m } symbols. 

This technique consists in estimating the carrier 
frequency f 0 of the transmitter or the complex 
z 0 =exp ( j2nf 0 T e ) in order thereafter to deduce the 
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symbols {a m } from the symbols {b m }, taking, from 
equation (3) : 



a m =b m exp(-j2nf 0 mIT e ) 



-ml 



(30) 



5 This step is applied after step 1.4 of reordering the 
symbols and the channel vectors. From equations (3), 
(17), (7) and (8), the following channel vectors are 
used: 



10 



h f (ii) = 



Z ;' h(nIT e ) 
z 0 nUi h{nlT e + T e ) 



*o -+(| -'ty«n;+(/-ra. 



(31) 



for neQ 



where Q={ Ind (rp - L o, rp +L0](n)=l for a p such that 1 £ p ^ P} 



15 



Since it is known that Q = {ni<...<n Lc } , a grand vector b 
is obtained from the vectors h z (n), such that: 



w = 



M"i) 

M"*)j 



(32) 



20 



The search for f 0 consists in maximizing the following 
criterion : 



Carrier (f 0 )=|w H c(exp(j27tf 0 T e ))| 2 



(33) 



where c(zo)= 



c(rt„z 0 ) 
c(w 2 ,z 0 ) 

c(n K ,z 0 \ 



and where c(n,Zo) = 



nt+\ 

Z 0 



nl +(/-!) 
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The steps of the process suitable for the case of a 
transmitter with a non-zero frequency are the 
following: 

5 Step III.a.1: Step 1.1 to step 1.4 described above in 

order to obtain the symbol trains associated with 

the channel vectors h z (kj). 

Step III. a. 2: Construction of the vector w of equation 

(32) from the fi z (kj) . 
10 Step III. a. 3: Maximization of the carrier (fo) 

criterion of equation (33) in order to obtain f 0 . 

Step III. a. 4: Application of equation (30) in order to 

deduce the symbols {a m } from the symbols {b m } . 

Step III. a. 5: Step 1.5 to step 1.7 described above. 
15 In the case of a transmitter with non-zero 

frequency and for decorrelated multipaths, the steps 

are the following: 

Step IZI.b.l: Step II. a. 1 to Step II. a. 4 described 
above in order to obtain the vector z(t) of equation 
20 (22). 

Step III.b.2: Application of the ICA methods [4], [5], 
[15] and [19] in order to estimate L c symbol trains 

{b m .j} associated with the channel vectors h Z/ j. 

Step III.b.3: Reordering of the symbol trains {b m%J } 
25 and of the channel vectors h z ,j applying steps A.l, A. 2 

and A. 3 in order to obtain the symbol trains {b m _ k } 

associated with the channel vectors h z ,j = fi z (kj). 

Step III.b.4: Construction of the vector w of equation 

(32) from the h z (kj) . 
30 Step III.b.5: Maximisation of the carrier (fo) 

criterion of equation (33) in order to obtain f 0 . 

Step III.b.6: Application of equation (30) in order to 

deduce the symbols {a m } from the symbols {b m } . 

Step III.b.7: Choice of the symbol train associated 
35 with the higher-modulus vector h z (i), namely, {a m _i} . 

Step III.b.8: Step II. a. 7 to step II. a. 9 described 

above . 



In the case of a transmitter with non-zero frequency 
and for correlated multipaths, the steps are for 
example the following: 

Step IIZ.o.l: Step Il.b.l to step II. b. 6 No. 2.2 in 
order to obtain the vector z(t) of equation (29). 
Step III.c.2: Application of ICA methods [4], [5], [15] 
and [19] in order to estimate the L c symbol trains 

{b m j} associated with the channel vectors fi 2/ j- 

Step III.c.3: Reordering of the symbol trains i^mj} 

and of the channel vectors h 2j applying steps A.l, A. 2 

and A. 3 so as to obtain the symbol trains {b m _ kj } 

associated with the channel vectors h zj = h z (kj) . 

Step III.c.4: Construction of the vector w of equation 
(32) from the h,(* y ) . 

Step III.c.5: Maximization of the carrier criterion 
(f 0 ) of equation (33) in order to obtain f 0 . 
Step III.c.6: Application of equation (30) for deducing 
the symbols {a m } from the symbols {b m } . 

Step III.c.7: Choice among the symbol trains of that 
one which is associated with the higher-modulus vector 
h z (i), namely {a m _i} . 

Step III.c.8 : Step II. b. 8 to step II.b.10 described 
above . 
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